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Abstract. We extend a Gaussian model for the internal electrical potential of a 
two-dimensional Coulomb gas by a non-Gaussian measure term, which singles out 
the physically relevant configurations of the potential. The resulting Hamiltonian, 
expressed as a functional of the internal potential, has a surprising large-scale 
limit: The additional term simply counts the number of maxima and minima of 
the potential. The model allows for a transparent derivation of the divergence 
of the correlation length upon lowering the temperature down to the Kosterlitz- 
Thouless transition point. 



PACS numbers: 6460, 6130J, 6470 

1. Introduction 

The two-dimensional Coulomb gas has found an abundant number of applications, 
ranging from the melting of two-dimensional crystals jj], || , vortices in superfluid films 
H and thin superconductors j| and arrays of Josephson contacts || to topological 
defects of thin liquid-crystal films ||. Common to these systems is the existence 
of two types of topological, point-like defects (positive and negative 'charges'), which 
interact according to the two-dimensional Coulomb law, i.e. the potential energy grows 
logarithmically with distance. At low temperatures positive and negative defects are 
bound together and form neutral dipoles, which start to unbind at a certain critical 
temperature - the Kosterlitz-Thouless transition point [7], g|. In fact, an energy of the 
order of AE ~ log A (A is the system size) is needed to break up two bound charges. 
On the other hand, the unbinding process implies a gain of entropy, which is also of 
the order of log A. Above the critical temperature entropy wins and a plasma-like 
gas of free charges is formed. The high-temperature phase is characterized by the 
existence of a screening length: The Coulomb interaction is shielded beyond the so 
called Debye-Hiickel length due to the formation of a cloud of opposite charges around 
a (test-) charge. This phase is well described by a Gaussian model for the continuous 
charge density p (Debye-Hiickel model) with the Hamiltonian 

K/T=^- J d 2 idV p{x)G(x-x')p{x')+ j ^ x Yz> (L1) 

where G r ~ — (27r) _1 log(r/^)+const (£ is a microscopic lengthscale) is the Green 
function of the two-dimensional Laplacian (— V 2 G(r) = S(r)), Ka is the charge-charge 
coupling strength and the mass z measures the distance to the transition temperature 
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z ~ T — T c and controls the density of charge pairs. The charge-charge correlation 
function reads (in Fourier space) 

C<P) « f 1 - -5X^-1 (L2) 
V p 2 +K A zJ 

where the screening length (Debye-Hiickel length £dh) is given by Idh — {Kaz)^ 1 ^ 2 . 
The Gaussian functional ([Ll]) applies to temperatures well above the transition 
temperature - it does not predict the essential singularity of the correlation length 
£ ~ exp(fe/ \JT — T c ) closer to the transition point. 

The aim of this paper is to propose a novel model Hamiltonian for the high 
temperature phase of the two-dimensional Coulomb gas, describing the fluctuations 
of the internal electrical potential. Contrary to the sine-Gordon theory, which has an 
unphysical auxiliary variable as its fundamental field, our model is based on a physical 
quantity (the potential) and, therefore, leads to a more intuitive understanding of the 
Kosterlitz-Thouless transition. 

The resulting continuum model displays the correct divergence of the correlation 
length in the vicinity of the critical point, as will be demonstrated within a variational 
calculation and, in addition, within renormalized perturbation theory. 



2. The model 



We start with the (Gaussian) Debye-Hiickel model (1.1) and introduce the internal 
electrical potential cj> via — V 2 = p. The partition function for the model reads 

D[cb] W[0]exp (-^ J A 2 x{V4>f - | J d 2 x(V 2 0) 2 ) (2.1) 

with the (trivial) measure W[(f>] = 1. Can we do better and find an appropriate 
measure W[(f>] for the internal potential, which singles out the relevant configurations 
and, therefore, go beyond the Gaussian case? Here is a suitable candidate: 



W[<f>] = exp (~ J d 2 x(V 2 0) 2 ^(V0) 



exp(-Hi/T) (2.2) 



with the Gaussian ('delta function') 6\(E) = (2n\y 1 exp(-£' 2 /(2A)). 

First, we investigate the mathematics of this ansatz. For A — > oo the measure 
W[4>] is Gaussian and simply renormalizes the mass z. For small A, on the other hand, 
the function 5\(\7(j)) picks up zeros of V0, i.e. local maxima, minima and saddle 
points of the potential <f>. We expand <j) in the vicinity of an extremal point x (where 
V4>(xq) = 0) up to second order in £ = x — Xq: 

0(x o +£) « 4>(xo) + {l/2)a l3 ^ 3 i,j = 1,2 

(summation over double indices is assumed) and find for the contribution of this 
extremum (for A — > 0) 



Wi/r = | J d 2 x(v 2 0) 2 ^(V0) 

- f d 2 f (a ) 2 1 cxd ( a ^ a ' k ^^ k 

4 J 4 [an) 2nX P I 2A 

2 



U 



/\vim\ (2.3) 
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where 77^2 are the eigenvalues of the second derivative didjip = a™ right at the extremal 
point. Maxima/minima yield a contribution TL\/T > u, since r\\ ■ r\2 > 0. For a 
symmetric maximum/minimum rji — 772 and consequently TL\/T — u; on the other 
hand, a 'symmetric' saddle point yields zero, since 771 = — r]2- Therefore, the nonlinear 
term acts as a kind of chemical potential for the number of extremal points [II], p"l||f|: 

7i\/T > u x # (local maxima and minima of 4>). (2-4) 

To summarize this point, the nontrivial measure W[</>] favours potentials with only 
few maxima and minima, which are in fact the physically relevant ones. Consider 
a typical charge configuration of a Coulomb gas system right above the transition 
temperature. Apart from a large number of dipols, which give rise to an effective 
dielectric constant, we find few free charges qi = ±1 at positions tv The corresponding 
potential <j>{r) — J2i Q<G(r — r,) has (few) maxima and minima at the locations of 
the charges ri (if we use a short-distance regularized Green function G, otherwise the 
potential <fi would diverge at the locations r^) and symmetric saddle points elsewhere, 
since A<fi(r) = for r/r,. 

We rescale the potential y/JIcj) — * <fi and the constant /iA — > A, define the mass 
r = Ka/ and obtain as the final model 

H/T = J d 2 .*Q(V 2 0) 2 + ^(V0) 2 + J(V 2 0) 2 <5a(V0)) . (2.5) 

The coupling constants u, A are dimensionless, whereas r is a relevant coupling with 
dimension r ~ length~ 2 -it is a measure for the deviation from the critical point 
t~T-T c . 



3. Variational approach 



We calculate an upper bound of the free energy T of our model (2.5) with the help of 
T=- log J D[4>] exp(-W/T) < T v + (H - H v ) v (3.1) 

and the ansatz 

H v /T = I J d 2 x ((V 2 ^) 2 + c(V^) 2 ) (3.2) 

where A and u> are variational parameters, T v is the free energy with respect to 7i v 
and (. . .}„ denotes the corresponding average. In fact, the optimal Gaussian fit has 
the form given above — a more general ansatz is not necessary and would reduce to 
an expression such as (|Q|). To evaluate ((V 2 0(r)) 2 5\{V (j>{r))) v we note that V^(r) 
and V 2 </)(r) are uncorrelated Gaussian variables and, therefore (we drop the argument 
r from now on), 

,oo ^ , o „ r o cxp(-.E 2 /((V0) 2 ) -E 2 /{2\)) 

« V2 * )2 >- (3.3) 



ir((V0)2) 1! + 2»A 



The expectation value ((V0) 2 ) diverges for uj and, therefore, the width A of the 
Gaussian S\ becomes irrelevant, justifiying the introduction of a sharp delta-function 

I The number of saddles is equal to the number of maxima and minima for a fiat geometry. 
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(A — > 0) in the model ( |2.5| ). Up to the constant (H v ) v , we obtain the upper bound of 
the free energy per area / 

A 2 A 2 2 A 2 

8nf<[ dslog(A(s + u,))+ / ds 3 +r f ds \ 

Jo Jo As(s + uj) J As(s + uj) 

f A2 s 2 ( f A " s \ ~* 
+2u / ds— / ds— (3.4) 



o 



As(s + uj) \J As(s + uj) 
where the first term represents T v , the second and third term are the expectation 



value of the Gaussian part of (2.5) and the last term represents (Hi) v - A is the upper 
cutoff momentum and s = p 2 , irds = irpdp = d 2 p using polar coordinates. Next, we set 
A = 1 (equivalently we can introduce dimensionless couplings) , eliminate A (variation 
of the bound / with respect to A yields A = 1 + (r — uj) log(l + 1/uj)) and arrive at 

8tt/ < log (1 + (r - uj) log(l + l/w)) + log(l + uj) + uj log(l + 1/uj) 

+2u( - . 1 - . - -u ) . (3.5) 
Vlog(l + l/w) J 

It can be tested afterwards that uj and (r — uj) log(l + 1/uj) become small enough in 



the critical region to approximate log(l + (. . .)) «(...) in the first two terms of (3.5). 
The expression is especially simple for u = 1/2 

8 7 r / <rlog(l + lM + log(i ; iM (3.6) 

yielding the minimum log (1 + 1/uj) = 1/ \fr or 

wwexp^--^ (3.7) 

where we have used 1 + uj^ 1 w uj" 1 for r — > 0. ^ = uj^ 1 / 2 is the correlation length of 
the best Gaussian fit, hence an estimate for the correlation length as a function of r 

e~exp(^) (3.8) 

which is the celebrated essential singularity of the correlation length in the vicinity of 
the Kosterlitz-Thouless transition. For general u we obtain uj = exp(— ^/2uJt). 

4. Renormalized perturbation theory 



We have also studied the model Hamiltonian (2.1) within a simple (renormalized) 
perturbation expansion. Using the formalism presented in |], [Io| we have 
calculated the two-point vertex function ^(p), which is the reciprocal Fourier 
transform of the Green function (<f>(x)<f>(y)) , to lowest order in the coupling constant u. 
Details of the derivation can be found in the appendix. The vertex function displays a 
highly unconventional large-scale behaviour if compared to usual critical field theory. 
The resulting effective couplings are in fact finite in the limit r — > and fixed cutoff. 
Since an instability shows up as a mere artefact of the expansion, we have to introduce 
a (finite) shift of the mass r — > r + u At, where we treat u At J d 2 x(W(f>) 2 /2 along 
with the interaction term as a perturbation. For fixed spatial dimension d = 2, fixed 
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cutoff momentum A = 1 , a sharp delta function A — > and up to 0(w 2 ), the two-point 
vertex function reads (see appendix) 

r 2 (p) = A 0{{ (p 4 + T oaP 2 ) (4.1) 

with the effective couplings 

2u 

r off =t + uAt ^ (4.2) 

(log(l + 1/r)) 2 

27/ 

A ef( = 1 + — --. (4.3) 

log(l + l/r) V ' 

As a consequence of the 5(V0)— term and in contrast to ordinary field theories, the 
log— terms show up in the denominator and, consequently, the effective couplings 
remain finite in the critical limit t — > 0. However, at a particular r, provided 
At = 0, the effective mass becomes negative, signalling the breakdown of the naive 
perturbation theory. Even worse, the correction to the mass, divided by r itself, tends 
to infinity in the limit r — > . To absorb this divergence, we renormalize the mass by 
setting 

At = y . (4.4) 

(log(l + l/r)) 2 

Now, the effective mass reads r off = r and the deviation from the critical point is 
To = t + u At ps tt At for small r. We arrive at 

r = o (4.5) 

(log(l + l/r oi ,)) 2 

or r cff ~ exp ^— yj2u/Toj in agreement with the variational calculation. 



5. Conclusions 



We have proposed an alternative model for the high-temperature phase of the two- 
dimensional Coulomb gas. ft describes the fluctuations of the internal electrical 
potential with the help of a measure term, which favours potentials corresponding 
to configurations of only few charges. The range of validity of this model apparently 
extends down to the critical point — the model yields the correct singular behaviour of 
the correlation length while approaching the Kosterlitz-Thouless point from above, as 
shown within a simple (renormalizcd) perturbation expansion. In addition, a Gaussian 
variational approximation scheme turned out to be successful, in contrast to the sine- 
Gordon theory |l2[ , where an analogous Gaussian variational ansatz fails and yields a 
wrong singularity for the correlation length Jl3| . 

Several questions could not be addressed, for example how to extract the singular 
behaviour of the free energy or the universal critical superfluid density. In our opinion, 
the model deserves future investigations — thanks to its fascinating and unconventional 
properties. 
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Appendix A. Calculation of the two-point vertex function 

For fixed spatial dimension d = 2, A — » and up to 0(it 2 ), the two-point correlation 
function reads 

(0(*)0(3/)) = (0(*)0(y)) o (A.l) 
- /" d 2 Z ((«j/)) (2ArW( 2 )) 2 + (V 2 ^)) 2 *(W(*))) 

V(*)0(y) (2 At (V0(z)) 2 + (V 2 0(z)) 2 S (V0(z)))) q ) + 0( U 2 ) 

where (. . .)o denotes the (Gaussian) average with respect to the quadratic part 
Ho = (1/2) J d 2 x((V 2 0) 2 + r(V0) 2 ) of the Hamiltonian ([D|). As explained in the 
main text, we have introduced a mass shift r — > t+mAt and treat it Ar J d 2 x(V0) 2 /2 
as a counterterm (additional perturbation). The bare Green function reads 



4 



1 

p2 (jp. _|_ r *J 



The counterterm adds 



Go(p) = / d 2 a;(0( : r)0(O)) o exp(-ip.a ; ) = 357-^^. (A.2 

G c (p) = -u At Go (pfp 2 = -uAt 1 . (A.3) 

More involved are the contributions from the measure. We define the Gaussian 
variables 0i = <j>(x),<f>2 = <p(y),p = V 2 0(z),i£ = V0(z), the five-component vector 
M = (0i, 02, p, -E), and write down the distribution of M with the help of a Gaussian 
transformation jy] 



P(M) = J d 5 Mexp (-^ djMiMj + iM ■ Afj (A 



4) 



where = (MjJVfA Q denotes the correlation matrix. Explicitly, we have (we drop 
the subscript from now on) 



2 d / 5d 2 .Ed0id02d j od 2 J E 



(*) = (0(x)0(y) (VVW) *(V0(*)), 
1 

x exp (-1 ((0? + 2 ) <0 2 ) + ? <(V 2 0) 2 ) + £ 2 <(V0) 2 ) /2 

+ 20i0 2 (0(x)0(y)) + 2p0 x <0(:z)V 2 0(2)) + 2p0 2 <0(y)V 2 0(z)) 
+ 2^ • (0(x)V0(z)) + 2£0 2 ■ <0(y)V0(z))" 



i 0i0i + 0202 + p> + E ■ E)\ fafcp^E) (A.5) 



where the correlations (V 2 0(2)V0(;z)) vanish by symmetry. Next, we perform the 
trivial E integration and the E integration and obtain 

W = (2^n\(V^) /d<M0 2 dpd0 ld 2 dp 



x exp 



■5 ((#+#) <^ 2 )+P 2 ((V 2 0) 2 ) 
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+ 20 x 2 (0(a;)0(y)) + 2p0 x (<j>( X )V 2 ^(z)) + 2p0 2 (0(y)V 2 0(z)) 
- — (0i (0(*)V0(z)) + 2 (0(y)V#*)>^ 

+ i (0i0i + 0202 + ppy\ 4>i4>2P 2 - (A. 6) 

The variables 0i, 02, p are still Gaussian variables, with, however, modified correlations 
which can be found from (A. 6). Therefore, using Wick's theorem 

2 

tt((V0) 2 



(*) = r7m7T2T <0(*)V 2 0(z)) <0(y)V 2 0(*)) 



( {mm > - msm^m^) . (A . 7) 



7T<(V0) 2 ) ^ <(V0) 2 

One of these terms is cancelled by (see chapter ||) 

(0(^)0^)) <(V 2 0) 2 <5(V0)) = (0(^)0(y)) (A.8) 
yielding the following contribution of the measure term: 

Gl(p) = ~2^((V0) 2 ) ((p 2 +t) 2 ~ ((V0)^ 2 )p 2 (p 2 Vr) 2 ) ' (A ' 9) 
The vertex function reads 

r 2 (p) = (Go + G C + Gi)- 1 = p 2 ( P 2 + t)+u Arp 2 
/ 4 ((V 2 0) 2 ) 2 \ n( a 



2vr((V0) 2 ) ^ ((V0) 5 
= A cff (p 4 + r otf p 2 ) (A. 10) 

with the effective coupling constants (up to order 0(u 2 ) and cutoff momentum A = 1) 

^ = l + -/(i/ d a p^ 7 ^=l + w ^ 1/ ,, (A.ll) 
ru ^<(V 2 0) 2 ) 



,2vr 4, <A " ' p 2 + r; * ' log(l + l/r) 

and 



T«tt = t + u At 



2tt((V0) 2 ) 2vr((V0) 2 ) 



2u 

= t + uAt ^. (A.12) 

(log(l + 1/r)) 2 
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